Adsorption of titanium and titanium dioxide on graphene: n and p-type doping 
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Ah initio calculations within the density-functional theory formalism are performed to investigate 
the ground state, electric charge doping, and electronic properties of titanium and titanium dioxide 
monolayers adsorbed on a graphene surface. A new ground state structrure of Ti monolayer adsorbed 
on graphene is reported which is shown to be stable up to T = 500 K. Effects due to lower and higher 
Ti adatoms coverage are studied. We find that the adsorbed Ti provides a strong n-type doping 
which supports recent experimental observations. On the other hand, Ti02 can induce both p- and 
n-type doping in the carbon monolayer depending on whether oxygen or titanium atoms are closest 
to the substrate. We identify the structures which are responsible for the experimentally observed 
autocompensation mechanism that leads to the reversion of adsorbate effects after oxidation of the 
adsorbed Ti. 

PACS numbers: 81.05.ue, 61.48. Gh, 68.43.-h, 68.43.Bc, 68.43.Fg 



I. INTRODUCTION 

The interaction of graphene with metal adsorbates is 
one of the important research topics since the first suc- 
cessful experimental exploration of graphene. [1] Virtu- 
ally any device assembly incorporating graphene will nec- 
essarily include graphene-metal contacts. 

Titanium has been widely used for metal contacts on 
graphene with good adsorption properties. [2] Moreover, 
the interaction between titanium and graphitic surfaces 
and nanotubes has attracted a lot of attention. [3 Ti 
atoms are able to act as adsorption centers for molecu- 
lar binding. [4 Furthermore, Ti coated nanotubes and 
graphene have also been proposed for high-capacity hy- 
drogen storage solution. [5, 6 

However, a recent DFT study, ^ as well as a very re- 
cent experimental observation, [8] showed strong struc- 
tural changes of the adsorbed Ti layer on graphene due 
to the interaction with oxygen. Experimentally, it was 
shown that the adsorbed titanium on graphene leads to 
substantial n-type doping and a reduction of graphene's 
mobility. The subsequent exposure of the samples to oxy- 
gen has restored graphene's gate dependent conductivity 
to almost intrinsic values, effectively cancelling any pre- 
vious metal-graphene interaction. [8 

Motivated by these experimental observations we per- 
form ah initio calculations to investigate titanium and 
titanium dioxide monolayers adsorbed on graphene. For 
the identified ground state structures, the electronic band 
structure is calculated and a charge population analysis 
is performed to investigate the character of the charge 
carrier doping in graphene. We find a strong n-type dop- 
ing of graphene induced by adsorbed Ti and identify the 
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highest possible monolayer coverage with stoichiometry 
TiaCg . When Ti02 monolayer crystals are aligned on 
top of the graphene surface we observe p- or n-type dop- 
ing of graphene, depending on the nature of the atoms 
that are exposed closest to graphene with oxygen acting 
as an acceptor and titanium as a donor. 

This paper is organized as follows. First, we describe 
the computational details of our ah initio calculations, 
followed by an investigation of the ground state structure 
of a titanium monolayer on graphene. Next the charac- 
ter of the charge transfer between adsorbed Ti and the 
graphene substrate is examined. Further model titanium 
dioxide monolayer structures are investigated. To con- 
clude, we report titanium dioxide structure dependent 
doping of graphene and link this observation with the ex- 
perimentally observed recovery of the graphene gate de- 
pendent conductivity after the oxidation of the adsorbed 
Ti atoms. 



II. CALCULATIONS 

All our calculations were performed within the den- 
sity functional theory (DFT) formalism as implemented 
in the VASP package with usage of the local spin den- 
sity approximation (LSDA) for the exchange-correlation 
functional. We made use of the projector augmented 
wave method [9] and a plane-wave basis set with an en- 
ergy cutoff of 500 eV. The complete relaxation of atomic 
positions and supercell size was performed up to the level 
that the forces are smaller than 0.01 eV A~^. 

For all the used supercells, a sampling of the Bril- 
louin zone was done with the equivalent of a 28 x 28 x 1 
Monkhorst-Pack pTO^ /c-point grid for the graphene unit 
cell (containing two carbon atoms). Spin polarization 
was included in the calculations. 

Because periodic boundary conditions were applied in 
all three dimensions, the height of the supercell was set 
to 20 A, to include enough vacuum to minimize the 
interaction between adjacent layers, and dipole correc- 
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tions were used. All reported quantitative results of the 
charge transfer were calculated with the iterative Hirsh- 
feld charge population analysis. [11] 



III. RESULTS AND DISCUSSION 

The binding energies are calculated with respect 
to intrisic monolayer graphene and the isolated 
atom/molecule of which the binding energy is reported. 
The formation energy for titanium covered graphene is 
calculated in relation to intrinsic monolayer graphene and 
the energy the of Ti atoms in a free standing monolayer 
titanium crystal with hexagonal symmetry. 



A. Properties of titanium monolayer on graphene 

Properties of single Ti atom adsorption on a graphene 
surface were already investigated and the most preferable 
adsorption site for a Ti atom was found to be the hollow 
site with a binding energy of -1.58 eV per Ti atom. The 
top and bridge adsorption site were found to be energet- 
ically less favorable by 0.62 eV and result in metastable 
states. [T2I 




FIG. 1: Structure g2x2-3Ti: Ti monolayer on a graphene 
sheet, 3 Ti atoms (light) per 8 C atoms (dark, displayed as a 
mesh in schematic picture on the left). 

Also high density coverage of graphene with titanium 
has been previously investigated and some ground state 
structures were proposed which have stoichiometry TiCg 
and TiC2 with Ti atoms placed on hollow sites creat- 
ing monolayer metalic coverage. [I3l |14] However, we re- 
port here a different ground state structure. As written 
above, a single Ti atom placed on top of a free standing 
graphene surface is strongly bonded at the hollow site 
with a diffusion barrier high enough to prevent its mo- 
tion on the surface even at room temperature. This fact 
suggests that the most preferable high coverage structure 
for titanium on graphene will be the one in which above 
every single hollow site there is one Ti atom (i.e., one 
Ti atom per two C atoms). Even the lattice constants 
compare favourably for this structure with a mismatch 
of only 5 % (see Table 

We have calculated binding energies for a variety of 
structures with stoichiometry TiC2 and we also investi- 
gated the structure g2x2-3Ti (displayed in Fig. [T]) with 



a lower amount of Ti atoms per C atom (three Ti atoms 
per eight C atoms) and a larger lattice mismatch of -9 
% (Ti atoms are spaced further apart from each other in 
comparison to the titanium monolayer crystal). The cal- 
culated formation energy of -0.93 eV per Ti for g2x2-3Ti 
system was found to be lowest among all the investigated 
structures. The lowest formation energy for the TiC2 sys- 
tem was found to be -0.76 eV per Ti, for the structure 
consisting of adatoms placed above hollow sites with ev- 
ery fourth Ti atom displaced further from the graphene 
substrate (see Fig. [2|d) initial state). Substantial forma- 
tion energy advantage makes the g2x2-3Ti system virtu- 
ally the only plausible structure for titanium atoms on 
graphene at large coverage densities. 

We have been surprised by this finding that the struc- 
ture g2x2-3Ti has the lowest binding energy among stud- 
ied cases. This has motivated us to perform DFT based 
molecular dynamics simulations of systems with different 
Ti coverage densities to identify the ground state configu- 
rations that may have been overlooked and/or to support 
the aforementioned findings. 

The simulated systems consist of a 4x4 graphene su- 
percell with 2, 4, 8, and 16 Ti adatoms, as displayed in 
Figs. |2|a-d). At the temperature of 500 K the different 
systems were allowed to evolve to more favourable states 
within a time window exceeding 1 ps. 

The system with the highest concentration coverage 
(TiC2), with Ti atoms initially placed above every hol- 
low site, does not remain a monolayer crystal. However, 
something interesting happens: the Ti atoms rearrange 
themselves so that they create one monolayer on the 
graphene surface with the same structure as g2x2-3Ti 
while the redundant Ti atoms create a second layer. 

The systems with lower concentration of Ti atoms per 
C atoms (TiCg and TiC4, Figs, ^h) andj^Jc), respec- 
tively) have also undergone structural changes. The sep- 
aration of adsorbed Ti atoms was too small to prevent 
atoms from mutual interaction. While adatoms remained 
adsorbed as a monolayer we have observed the creation 
of clusters and local rearangements. Those structural 
changes occured despite the presence of the diffusion bar- 
rier for adatoms adsorbed on the hollow site (clearly vis- 
ible for TiCg, Fig. [2|b)). Again the comparison with the 
structure g2x2-3Ti reveals remarkable similarities in local 
arrangements of Ti adatoms, which is more pronounced 
in the structure TiC4 (Fig. [2|c)) with higher adatom 
coverage. 

The last of the systems investigated using MD calcula- 
tions is the one with the lowest Ti concentration - TiCie 
(Fig. ^8i)). We found that this case is the only one 
that maintained its structural properties during the MD 
simulations. As can be seen in Fig. |2|a) the adsorbed 
atoms are separated by two carbon hexagonal rings. The 
separation distance of 6.5 A is found to be sufficient to 
suppress the mutual adatom interaction observed in the 
previous cases. The Ti atoms were trapped at the hollow 
sites and do not pass across the diffusion barrier. 

The aforementioned results of the MD calculations 
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FIG. 2: Top and side view of different graphene structures 
with different amount of titanium coverage before (left) and 
after (right) 1 ps of relaxation at T = 500 K. The systems 
consist of (a) 2, (b) 4, (c) 8 and (d) 16 Ti adatoms (light 
symbols) placed on the graphene sheet in an unit cell with 32 
C atoms (dark symbols). 



can be summarized as follows: (i) Ti atoms placed on 
a graphene surface can overcome diffusion barriers (for 
the hollow site) by mutual close range interaction (up to 
±6.5A), (ii) the maximal titanium monolayer coverage 
of graphene has stoichiometry TiaCg and corresponds to 
the strucure g2x2-3Ti (Fig. [l]). 

Experimentally [8^ a strong charge transfer from the Ti 
layer to graphene was observed. Charge population anal- 
ysis performed on the investigated systems has revealed 
the same trend of strong to moderate charge transfer 
from the Ti atoms to the underlying graphene for low and 
high concentrations, respectively. However, the quantita- 
tive and even qualitative values for charge transfer were 



found to be inconsistent among the different methods 
we used (including the iterative Hirshfeld charge popu- 
lation analysis, modified iterative Hirshfeld method and 
the Bader and Voronoi cells charge population analysis). 

imiiMHi 

An alternative to obtain directly the charge transfer 
is by calculating the dipole moment after Ti adsorption. 
However, the value of the electric dipole moment for the 
g2x2-3Ti structure is 0.087 eA, suggesting a small charge 
transfer from graphene to Ti. This corresponds neither 
to the experimental observation of n-type doping, nor 
with the charge transfer analyses discussed above. As 
shown by the plane- averaged deformation electron den- 
sity in Fig. |3) the strong oscillations in the charge den- 
sity do not clearly indicate actual charge transfer nor a 
dipole that will identify such transfer. The explanation 
for opposite total electric dipole moment lies in the domi- 
nant polarisation of the carbon layer induced by the pos- 
itively charged chemisorbed Ti monolayer. The surplus 
of charge located on the C atoms is shifted towards the 
Ti atoms thus creating a dipole moment counteracting 
and exceeding the dipole moment based on the negative 
charge located on the carbon layer and the positively 
charged Ti atoms. We expect that the polarisation of 
the carbon layer to be the probable reason for the in- 
consistent results obtained by the aforementioned charge 
population analysis methods. 
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FIG. 3: Plane- averaged deformation electron density of g2x2- 
3Ti configuration obtained from the difference between the 
electron density of the g2x2-3Ti structure and the electron 
densities of the subsystems consisting of only carbon or tita- 
nium atoms. 



TABLE I: Lattice constants and workfunctions. 



structure 


$ (eV) 


a (A) 


graphene 


4.49 


2.447 


hexagonal monolayer Ti 


4.72 


2.566 


rectangular monolayer Ti02 


8.54 


3.141 


hexagonal monolayer Ti02 


8.51 


2.945 



Since the direct charge transfer calculations proved to 
be insufficient in determinig the type of doping we also 
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FIG. 4: Band structure of the g2x2-3Ti configuration. Tiie 
energies are relative to the Fermi level (Ef = 0). The Dirac 
point lies -2.61 eV below the Fermi level showing strong n- 
type doping. The amount of carbon p^ character is indicated 
by the blackness of the bands. 



performed electronic band structure calculations and cal- 
culated the workfunction of the pristine and hybrid ma- 
terials. As can be seen from the electronic band structure 
of the g2x2-3Ti configuration (see Fig. [4|, strong n-type 
doping occurs to graphene. The Dirac cone, still present 
in the electronic band dispersion, is shifted by -2.61 eV 
below the Fermi level. The calculation of the workfunc- 
tion reveals the same type of doping, albeit in a partially 
counterintuitive way. As shown in Table [T] the work func- 
tion of free-standing graphene is 4.49 eV, while the work 
function of an isolated hexagonal titanium monolayer is 
4.72 eV. From a simple comparison of these values there 
is no precondition for charge transfer to occur from Ti 
to the graphene substrate. However, for closely sepa- 
rated metal-graphene structures (2.1 A for g2x2-3Ti) the 
chemical interaction prevails and is responsible for n-type 
doping even when the workfunction of the metal is larger 
than the graphene work function. [19] In the final sys- 
tem the computed workfunctions on graphene (3.97 eV) 
and titanium side (4.72 eV) indicate a charge transfer to 
graphene. The workfunction difference increases as com- 
pared to the pristine values (0.81 eV compared to 0.23 
eV). 



B. Properties of Titanium dioxide monolayer on 
graphene 

Let us now consider the case of Ti02 adsorption on 
graphene. As a model approach for studying the in- 
teraction of titanium dioxide on graphene surfaces we 
have chosen monolayer Ti02 films derived from anatase, 
(101) with hexagonal and (001) with rectangular symme- 
try, which were proven to be stable. [20] From the broad 
range of possible configurations the structure g2x2-hex- 
Ti02 with lattice mismatch of about 4 % (displayed in 
Fig. [sjc)) was found to be the most stable. Ti atoms are 
aranged in a monolayer triangular lattice, similar to the 



one found in the previously discussed structure g2x2-3Ti, 
and oxygen atoms are placed between every three neigh- 
bouring Ti atoms. The Ti02 monolayer is buckled due 
to the alternating displacement of oxygen atoms above 
or below the plane of Ti atoms. The other displayed 
model configurations in Figs, [sja) andjsjb), while being 
metastable, will be discussed below as examples which 
lead to different doping properties. 



a) 



o.., 









x" / 



O' 





o ^ 




FIG. 5: Model structures for a Ti02 monolayer on a 
graphene sheet. The systems (a) g2x2- aligned hex-Ti02 
and (c) g2x2-hex-Ti02 contain 3 Ti02 molecules per 8 C 
atoms, whereas (b) the g-sq8-Ti02 configuration has 2 Ti02 
molecules per 8 C atoms. The carbon atoms are given by gray 
colour (a mesh in schematic picture on the left), Ti atoms are 
in light gray and O atoms are shown in red (dark, small circles 
in schematic picture on the left). 

The charge population analysis performed on the 
Ti02-graphene structures were found to be much more 
robust than for the titanium covered graphene. The Ti-0 
bonds exhibit an ionic character with Ti atoms gaining 
positive charge after loosing two electrons in favour of the 
O atoms. Besides the charge redistribution in the Ti02 
layer the g2x2-hex-Ti02 structure shows also p-type dop- 
ing of the graphene layer. The amount of charge trans- 
fered from graphene, calculated with iterative Hirshfeld 
charge population analysis, is 0.126 e per Ti02 unit. 

The character of the doping is also visible in the elec- 
tronic band structure shown in Fig. |6]with a clear hole 
doping of graphene and a shift of the Dirac cone above 
the Fermi level by 0.5 eV. The observed doping can be ex- 
plained from geometrical arguments. The bounded oxy- 
gen and titanium atoms can induce p- or n-type doping 
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FIG. 6: Band structure of g2x2-hex-Ti02 configuration. Tiie 
energies are relative to the Fermi level (Ef = 0). The Dirac 
point lies 0.5 eV above the Fermi level showing p-type dop- 
ing. The amount of carbon p^ character is indicated by the 
blackness of the bands. 

when close to the graphene surface. The Ti02 monolayer 
in the g2x2-hex-Ti02 structure can be seen as consisting 
of three layers - a titanium layer sandwiched between two 
oxygen layers, thus exposing oxygen atoms to graphene. 
Similar charge transfers have been observed at the in- 
terface between graphene and rutile Ti02 (110), where 
hole doping in graphene and an accompanying shift of 
the Fermi level relative to the Dirac cone was observed. 

m 

The aforementioned reasoning proves to be valid also 
for the other systems we have investigated. The model 
structures in Figs. |5|a) andjsjb) differ by the amount of 
oxygen or titanium atoms they expose to the graphene 
surface. The charge population analysis as well as the 
electronic band structure provide consistent evidence of 
p- or n-type doping depending on whether oxygen or tita- 
nium atoms are closer to the carbon atoms, respectively. 
The shift of the Dirac cone with respect to the Fermi 
level is -1.69 eV for g2x2-aligned hex-Ti02 and eV for 
g-sq8-Ti02 structure, and the amount of charge, calcu- 
lated with iterative Hirshfeld method, transfered from 
graphene is -0.452 e and 0.026 e, respectively. 

Recently it was found experimentally that the 
graphene gate dependent conductivity can be recovered 
after oxidation of adsorbed Ti atoms. [HJ This recov- 
ery can be explained with the autocompensation of the 
titanium dioxide surface [22 exposed to graphene. We 
showed above that Ti02 structures can induce both p- 



or n-type doping in graphene thus opening the possibil- 
ity of mutual compensation of the doping effects. Al- 
though the groundstate Ti02 monolayer induces hole 
doping to graphene the performed calculations suggest 
that metastable states of Ti02 monolayers can be re- 
sponsible for the recovery of the graphene gate depen- 
dent conductivity. However, the geometrical nature of 
the autocompensation does not imply the actual pres- 
ence of the above considered monolayer Ti02 crystals on 
the graphene surface in the experiment. Existence of any 
Ti02 nanostructure with both p- or n-type doping sides 
exposed to the graphene surface will lead to similar ef- 
fects. 



IV. CONCLUSIONS 

We have demonstrated strong n-type doping of Ti cov- 
ered graphene. While the various charge population anal- 
ysis methods failed to provide reasonable quantitative re- 
sults for the shift of the Dirac point with respect to the 
Fermi level our investigation of the increased difference 
of the work functions {^m — ^g) provided a consistent 
picture of the character of the induced doping. Titanium 
is found to be strongly chemisorbed on graphene and the 
highest monolayer coverage is found to have stoichiome- 
try TiaCs . The studied Ti02 monolayer crystals aligned 
on top of graphene induce charge doping that depends on 
the nature of the atoms that are closest to graphene with 
oxygen acting as an acceptor and titanium as a donor. 
The ground state structure of Ti02 monolayer covered 
graphene exhibits p-type doping of graphene, neverthe- 
less the calculations suggest these metastable states of 
Ti02 monolayer can be responsible for the recovery of the 
graphene gate dependent conductivity after oxidation of 
adsorbed titanium, as recently found in experiment. [8]. 
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